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Single  edge-notched  bend  (SENB)  specimens  containing  shallow  cracks  («/ 
W<02)  are  commonly  employed  for  fracture  testing  of  femoc  materials  to  (he 
lower-transition  region  where  extensive  plasbaty  (but  no  agwhcant  ductile  crack 
growth)  precedes  unstable  fracture.  Critical  ./-value*  (jfc)  for  shallow  crack  specs- 
mens  are  significantly  larger  (factor  of  2-3)  than  the  ^^vmhtes  for  correepondmg 
deep  crack  specimens  at  identical  temperatures.  The  increase  of  fracture  tough¬ 
ness  arises  from  the  loci  of  constraint  that  occurs  «kta  the  groni  ptasnc  aooes  of 
bending  impinge  on  the  otherwise  autonomous  crack-bp  plane  tones  Conse¬ 
quently,  SENB  specimens  with  small  and  large  s/Wrnoo*  loaded  to  (he  same  lva¬ 
lue  have  markedly  different  crack-bp  stresses  under  largn-ccale  pUsoary.  De¬ 
tailed,  plane-strain  finite-element  analyses  and  a  local  mti  w  bused  criterion  for 
cleavage  fracture  are  combined  to  establish  yeomen  site  requirements  (deforma¬ 
tion  limits)  for  testing  in  the  transition  region  which  secure  a  single  parameter  (S) 
characterization  of  the  crack-tip  strees  field.  Moreover,  these  analyses  provide  the 
first  quantitative  framework  to  correlate  <4 -values  with  */W  ratio  once  the  defor¬ 
mation  limits  are  exceeded.  The  new  procedures  are  dtreedy  applicable  to  eetimest 
the  toughness  values  for  shallow  crack  yeomens  from  remits  of  teas  performed 
on  conventional  deep-crack  yea  mem.  The  increased  toughness  values  for  dial* 
low-crack  specimens  are  required  to  perform  realotic  smrmcnts  of  the  surface¬ 
breaking  defects  most  commonly  encountered  in  strumites-  The  correlation  pro¬ 
cedure  is  applied  to  an  extensive  set  of  test  results  for  an  A36  steel  and  is  shown  to 
successfully  correlate  the  a/IKratio  effects  on  the  measured  fracture  coughoem  val¬ 
ues  in  the  transition  range. 
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EXECUTIVE  SUMMARY 


This  study  provides  the  first  quantitative  step  toward  development  of  a  framework  iu 
address  the  issue  of  constraint  that  severely  hampers  the  aacsamem  of  shallow,  ua 
face-breaking  defects.  Realistic  assessments  require  that  crack-up  amdtuom  (u 
stress  and  strain)  in  simple  laboratory  test  specimens  closely  match  those  at  the  defect 
in  a  structural  component.  When  similar  crack -op  coodtuam  prevail,  critical  values  of 
macroscopic  toughness  (/)  measured  in  the  test  yeomen  are  directly  applicable  to  «* 
sess  the  criticality  of  the  structural  defect.  However,  large-scale  plastic  deformation 
and  effects  of  geometry  lead  to  complex  relaoonshtpt  between  J and  v  -nan -strum  «t 
the  crack-tip  that  are  not  currently  understood 

Fracture  toughness  values  obtained  from  extensive  testing  of  shallow-notch  bend 
specimens  have  been  shown  previously  to  correlate  qualitatively  with  the  fracture  o! 
specimens  containing  shallow,  surface  defects.  The  absence  of  a  quantitative  undet 
standing  for  this  behavior  motivated  the  current  study  The  effects  of  tn -plane  plasm 
flow  on  constraint  for  through-cracks  in  bend-bars  are  essentially  revolved  in  the 
present  work  through  careful  numerical  analyses  and  comsderatson  of  a  single  fract  vie 
mode  (stress-controlled).  Quantitative  guidelines  to  estimate  fracture  toughness  val¬ 
ues  for  any  alW ratio  given  test  results  at  one  a/W rapo.  generally  0.5.  are  an  tmmedi 
ate  benefit  of  these  numerical  computations.  Moreover,  the  procedure  developed 
here  to  quantify  in-plane  constraint  appears  equally  promising  to  quantify  thickness 
effects  for  through-cracks  in  bend  bars  which  leads  directly  to  the  surface  crack  con 
figuration  of  most  practical  interest. 
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ABSTRACT 


Single  edge-notched  bend  (SENB)  specimens  containing  shallow  cracks  (a! 
W<0.2)  are  commonly  employed  for  fracture  testing  of  ferritic  materials  in  the 
lower-transition  region  where  extensive  plasticity  (but  no  significant  ductile  crack 
growth)  precedes  unstable  fracture.  Critical  /-values  (/c)  for  shallow  crack  speci¬ 
mens  arc  significantly  larger  (factor  of  2-3)  than  the  /c-values  for  corresponding 
deep  crack  specimens -at  identical  temperatures.  The  increase  of  fracture  tough¬ 
ness  arises  from  the  loss  of  constraint  that  occurs  when  the  gross  plastic  zones  of 
bending  impinge  on  the  otherwise  autonomous  crack-tip  plastic  zones.  Conse¬ 
quently,  SENB  specimens  with  small  and  large a/W'ratios  loaded  to  the  same/-va- 
lue  have  markedly  different  crack-tip  stresses  under  large-scale  plasticity.  De¬ 
tailed,  plane-strain  finite-element  analyses  and  a  local  stress-based  criterion  for 
cleavage  fracture  are  combined  to  establish  specimen  size  requirements  (deforma¬ 
tion  limits)  for  testing  in  the  transition  region  which  assure  a  single  parameter  (/) 
characterization  of  the  crack-tip  stress  field.  Moreover,  these  analyses  provide  the 
first  quantitative  framework  to  correlate  /c-values  with  alW ratio  once  the  defor¬ 
mation  limits  are  exceeded.  The  new  procedures  are  directly  applicable  to  estimate 
the  toughness  values  for  shallow  crack  specimens  from  results  of  tests  performed 
on  conventional  deep-crack  specimens.  The  increased  toughness  values  for  shal¬ 
low-crack  specimens  are  required  to  perform  realistic  assessments  of  the  surface- 
breaking  defects  most  commonly  encountered  in  structures.  The  correlation  pro¬ 
cedure  is  applied  to  an  extensive  set  of  test  results  for  an  A36  steel  and  is  shown  to 
successfully  correlate  the  a /IV ratio  effects  on  the  measured  fracture  toughness  val¬ 
ues  in  the  transition  range. 
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INTRODUCTION 


Single  edge-notched  bend  (SENB)  specimens  containing  shallow,  through  cracks  (at 
W<  0.2)  are  commonly  employed  for  fracture  toughness  testing  in  the  ductile-to-brittle 
transition  region  of  ferritic  materials.  The  testing  of  particular  microstructures  in  weld¬ 
ments  frequently  dictates  the  use  of  shallow  crack  specimens  when  deeply  cracked  speci¬ 
mens  are  impossible  to  obtain  due  to  material  thickness  and  weld  orientation.  Shallow 
crack  specimens  also  provide  an  opportunity  to  better  match  the  crack-tip  conditions  in  a 
structural  component  containing,  for  example,  a  surface  flaw  loaded  primarily  in  tension. 
Fracture  toughness  values  obtained  from  deeply  notched  bend  specimens  characterize 
microstructures  and  constraint  conditions  rarely  experienced  by  in-service  structures, 
making  such  specimens  inappropriate  for  use  in  assessments  of  structural  integrity. 

Recent  experimental  studies  [  1-9]  have  demonstrated  a  significant  increase  in  the  ap¬ 
parent  elastic-plastic  fracture  toughness,  characterized  by  the  /-integral  [10],  for  speci¬ 
mens  with  shallow  cracks  when  tested  in  the  transition  region.  Figure  1  illustrates  the  dif¬ 
ferences  in  toughness-temperature  relationships  typical  of  a  mild  structural  steel  tested 
with  deep  ( alW  *  0.5)  and  shallow  (alW  *  0.15)  crack  specimens.  Load-displacement  re¬ 
cords  for  specimens  tested  throughout  the  transition  region  show  some  degree  of  non¬ 
linearity;  plastic  hinges  commonly  form  prior  to  failure  even  towards  the  lower-end  of  the 
transition  region.  At  upper-shelf  temperatures,  ductile  fracture  initiation  by  microvoid 
coalescence  and  shear  localization  is  generally  linked  to  the  achieving  of  a  critical  strain 
over  an  intensely  deformed  process  zone  at  the  blunted  crack  tip  [11].  Conversely,  at  or 
near  the  lower-shelf  temperatures,  brittle  (unstable)  fracture  initiation  by  cleavage  is  at¬ 
tributed  to  achieving  a  critical  tensile  stress  over  a  significant  microstructural  distance 
[12,13].  The  initiation  of  stress-controlled  cleavage  fracture  in  an  SENB  specimen  de¬ 
pends  upon  material  properties  (cleavage  fracture  stress,  flow  characteristics)  and  struc¬ 
tural  parameters  which  govern  crack-tip  conditions  ( alW  ratio,  specimen  thickness). 
Assessments  of  structural  integrity  generally  focus  on  the  cleavage  fracture  mechanism 
due  to  the  possibility  of  a  local  crack-tip  instability  triggering  a  catastrophic  structural 
failure. 

In  the  lower-transition  region,  unstable  fractures  with  cleavage  as  the  microstructural 
separation  process  occur  following  plastic  deformation  that  exceeds  A/c  validity  limits 
(ASTM  E3 99-83)  but  without  significant  ductile  tearing.  The  onset  of  stable,  ductile  tear¬ 
ing  represented  byJjc  (ASTM  E813-87)  defines  the  upper-transition  region,  although  un¬ 
stable  fractures  triggered  by  cleavage  are  sometimes  observed  after  ductile  tearing.  Criti¬ 
cal  /-values  for  cleavage  fracture  toughness  in  the  lower-transition  region  are  denoted  Jc 
to  differentiate  from  the  larger  values  (//c)  that  apply  to  the  onset  of  ductile  tearing.  As 
shown  in  Figure  1,  /c-values  for  short  crack  specimens  may  be  significantly  larger  than 
values  for  deep  crack  specimens  at  identical  temperatures  in  the  lower-transition  region. 
The  increased  ductility  of  shallow  crack  specimens  also  lowers  the  transition  temperature 
and  increases  the  ductile  initiation  toughness,  //c  [1,6-8]. 
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The  increase  of  /<  -values  in  shallow  notch  tpcomcm  develops  when  the  u» -plane  pu>- 
tic  deformation  of  the  gross  bending  impinges  on  the  local  crack -up  fields  which  itluri 
the  kinematic  constraint  against  funher  pUM*.  flow  The  yutmag  of  a  at*,  up  pUviu  runes 
with  those  at  the  tension  surface  toe  shallow  notches  vharpb  reduces  the  mr»  uuuualm 
(see  Fig.  1).  Deeply  notched  specimens  of  very  low  hardening  matrnah  also  edutui  «  k*» 
of  crack-tip  constrauu,  but  not  as  severe  as  foe  shallow  notches  Once  thr  global  and  local 
plastic  fields  interact,  the  crack -up  strains  and  stresses  no  longer  increase  m  proportion  n> 
one  another  according  to  the  asymptotic  singular  fields  derived  by  Hutchinson  1 1  a]  Hue 
and  Rosengren  (15)  (HRR)  which  have  amplitudes  uniquely  governed  by  the  single  pa 
rameter  J.  The  shallow  crack  specimens  edubsi  smaller  stresses  and  larger  strains  than 
those  represented  by  the  highly  constrained  HRR  fields,  thus  deep  and  shallow  notched 
specimens  loaded  to  the  same  applied  /-value  (that  produces  significant  pUuuatv  t  cshibii 
very  different  stress  fields  at  the  crack-up  which  leach  to  the  specimen  utr  dependence  ot 
critical  /c-values.  However.  J  retains  usefulness  as  a  m*c*o*cup»c  measure  of  fractutc 
toughness,  albert  geometry  and  material  dependent,  to  link  the  crack = up  stress  strain 
fields  with  imposed  loading  even  though  strict  coupling  with  the  HRR  fields  n«»  longer 
holds. 


McMeckmg  and  Parks  1 16).  Shih  and  German  J 17).  and  Stub  )  18)  performed  plant 
strain  finite -clement  analyses  for  deeply  notched  SEND  speamem  to  evublnh  the  levels 
of  plastic  deformation  beyond  which  the  single  parameter  HRR  fields  no  longer  dominate 
or  characterize  the  complete  crack -up  corsduiom  The  requirement  may  be  Mated  as 


b  2 


U<r  Jk 

Gfkm 


f*<r 


25  for  a/W  *  0  5 


n  > 


where  b  is  the  uncracked  ligament  length  fW-a)  and  o**,  n  the  material  Arm  micks.  dc 
fined  as  the  average  of  the  yield  and  tensile  strength  Generally,  this  'unit  is  satisfied  at 
deformations  well  into  the  large-scale  yielding  range  represented  by  the  formation  of  a 
plastic  hinge.  To  establish  the  value  for  n<, .  general  agreement  ( •  10%)  between  the  fim 
te-element  and  HRR  crack-tip  fields  was  sought  over  a  distance  of  3- 5  <5, .  where  <5,  dc 
notes  the  crack  tip  opening  displacement  (CTOD)  No  elicit  consideration  was  given  of 
a  fracture  separation  process  in  assigning  a  value  for  .  although  the  intended  use  was 
for  Ju  testing.  Parks  [19.20]  notes  the  inherent  difficulty  of  such  approaches  to  establish 
limits  assuring  /-dominance  because  the  strew  fields  change  gradually  from  the  /-dorm 
nated  ones  as  loading  is  increased,  assigning  values  for  iq,  thus  introduces  subjective 
judgements.  Nevertheless,  Eq  (1)  has  proven  adequate  to  guarantee  nearly  size-inde¬ 
pendent  Jic  values  for  initiation  of  ductile  tearing  and  is  incorporated  in  the  ASTM 
E8 13-87  test  procedure.  However,  this  requirement  is  now  believed  too  lenient  for  unsta¬ 
ble  fractures  resulting  from  stress-controlled  cleavage  because  even  a  slight  reduction  of 
crack-tip  stresses  below  the  HRR  values  can  produce  a  significant  elevation  of  the  critical 
/-value  (21-23). 
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The  present  work  focuses  on  quantifying  specimen  size  effects  for  unstable  fracture  by 
cleavage  when  plastic  deformation  exceeds  the  customary  Ku  validity  limits.  Very  detailed 
firute-element  analyses  along  the  lines  of  those  in  ( 17,18]  are  generated  for  both  shallow 
and  deep  notch  bend  specimens  for  a  wide  range  of  strain-hardening  behavior.  The  fini¬ 
te-element  models  fully  resolve  crack-tip  stress  fields  over  the  length  scale  of  3- 10<5,  for 
the  full  range  of  loading.  A  microscopic,  stress-based  criterion  for  cleavage  fracture  is 
adopted  to  establish  specimen  size  requirements  (deformation  limits)  which  assure 
unique  characterization  by  a  macroscopic  parameter  (J)  at  failure.  The  recommended  size 
requirements  are  significantly  less  stringent  than  those  of  E399-83  but  more  restrictive 
than  the  E8 13-87  requirements  of  Eq.  (1).  Moreover,  by  requiring  identical  values  of  the 
microscopic  criterion  at  failure,  the  finite-element  analyses  provide  a  framework  to  cor¬ 
relate  Jc -values  with  atW  ratio  once  the  deformation  limits  are  exceeded  (the  microme¬ 
chanism  of  tailure  is  assumed  to  remain  cleavage).  The  correlation  procedure  is  shown  to 
successfully  remove  the  geometry  dependence  of  fracture  toughness  values  in  the  lower- 
transition  range  for  an  A36  steel. 


STRESo  FIELDS  FOR  ANALYSIS  OF  CLEAVAQE  FRACTURE 


A  study  of  specimen  size  effects  on  cleavage  fracture  toughness  must  correlate  the  various 
crack-tip  stress  fields  over  the  potential  fracture  process  zone.  The  plane-strain  HRR 
fields  provide  the  starting  point  for  such  comparisons  and  are  given  by. 
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where  r  and  6  are  the  polar  coordinates  centered  on  the  crack  dp,  /„  is  an  integration  con¬ 
stant  (function  of  n  only),  andd^  and  i y  are  dimensionless  functions.  Tables  of  the  HRR 
functions  are  provided  in  [24],  The  material  follows  a  pure  power-law  constitutive  model 
which  in  uniaxial  form  is 


r/€o  «  a(o/Oof 


(3) 


where  a0  an  effective  yield  stress,  «o  *  o0/£  is  the  reference  (yield)  strain;  a  is  a  dimen¬ 
sionless  factor  and  n  is  the  strain  hardening  exponent.  Generalization  of  this  relation  to 
multiaxial  states  using  deformation  (/?)  plasticity  yields  strain  components  of  the  form 
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where  jy  is  the  stress  deviator  and  a,  is  the  Mises  equivalent  tensile  stress. 

Eqs.  (2-4)  provide  the  asymptotic,  small-strain  theory  results  for  an  ideally  sharp 
crack-tip  under  monotonically  increasing,  Mode  I  loading.  These  fields  are  not  applicable 
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within  the  finitely  deformed  blunting  zone  which  evolves  at  the  crack-dp.  The  crack-up 
opening  displacement,  <5, ,  provides  a  convenient  length  scale  to  describe  near-up  dis¬ 
tances,  such  as  the  blunung  zone  size.  Shih  (25)  defined  <5,  as  the  separation 
where  ±  45°  lines  emanadng  from  the  crack  dp  intercept  the  opening  faces.  Foi  the  HRR 
displacement  fields,  the  form  of  <5,  is 

6,  =  d„(a€0,n)J/a0  (5) 

where  the  coefficient  dn  is  weakly  dependent  on  aco  and  varies  from  0.3  for  small  n  to  0  8 
for  large  n. 

The  length  scale  over  which  the  single-term  fields  in  Eqs.  (2)  adequately  describe  the 
response  at  finite  distances  from  the  crack  dp  is  obtained  through  consideration  of  the 
boundary-layer  problem  referred  to  as  small-scale  yielding  (SSY).  Originally  proposed 
by  Rice  and  Tracey  [26],  McMeeking  (27]  and  recently  used  by  Parks  and  Wang  1 19].  the 
SSY  model  consists  of  an  annular  region  containing  a  sharp  crack  tip  which  is  subjected  to 
increasing  displacements  of  the  elastic  (Mode  I)  singular  field  on  the  outer  circular  bound¬ 
ary.  Steady  state  conditions  develop  in  which  strains  and  stresses  at  all  angles  scale  with 
r/(J/a0)  and  r/6, ,  as  do  the  HRR  fields,  and  persist  under  increasing  deformation  until 
the  plastic  zone  senses  the  circular  boundary.  Finite-element  analyses  for  the  SSY  prob¬ 
lem  indicate  diveigence  of  the  HRR  and  SSY  fields  for  r  >  2-3  6, ,  with  smaller  stresses 
and  larger  strains  predicted  by  the  SSY  model.  Thus,  Eqs.  (2)  may  be  considered  the  lead¬ 
ing  term  of  the  complete  near-tip  fields  provided  by  the  SSY  model. 

Finite  deformation  (blunting)  analyses  for  the  SSY  model  and  for  deeply  notched 
bend-bar  models  by  McMeeking  [27]  and  McMeeking  and  Parks  [16]  predict  the  location 
of  maximum  tensile  stress  on0  =0  at  approximately  2-46,  ahead  of  the  blunted  tip. 
Opening  stresses  nearer  to  the  crack  tip  are  relaxed  by  a  loss  of  triaxiality  due  to  the  free 
surface  of  the  blunted  tip.  The  fracture  process  zone  for  cleavage  thus  lies  beyond  the  fin¬ 
itely  deformed  zone  adjacent  to  the  blunted  tip  and  beyond  the  region  over  which  the 
asymptotic  HRR  field  is  adequate  to  describe  the  stress  magnitudes.  Moreover,  differ¬ 
ences  in  stresses  for  the  SSY  model  predicted  by  the  notch-blunting  analyses  and  by  the 
small-strain  analyses  (with  a  sharp  crack)  become  insignificant  for  r  >  2-3  6, .  These  ob¬ 
servations  suggest  that  detailed  finite-element  solutions  for  SENB  specimens  based  on 
small-strain  theory  are  valid  to  assess  fields  ahead  of  the  crack  tip  when  stress-controlled 
cleavage  is  the  fracture  separation  process  of  interest.  Finite-element  meshes  must  fully 
resolve  the  stress  fields  over  a  length-scale  of  2-10  6,  for  the  levels  of  loading  examined. 
To  assess  effects  of  specimen  size,  the  SENB  stresses  over  this  region  should  be  compared 
to  those  of  the  SSY  model  rather  than  HRR. 

With  the  basis  for  comparison  of  crack-tip  stress  fields  established,  a  simple  micro¬ 
scopic  model  for  cleavage  initiation  is  adopted.  Unstable  fracture  is  defined  to  occur  when 
the  tensile  stress  on  the  crack  plane  at  r = 46,  reaches  a  critical  value.  In  subsequent  sec¬ 
tions,  the  particular  choice  of  4,  5,  or  66,  is  shown  to  be  insignificant  as  the  stress  fields 
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scale  identically  with  macroscopic  loading  parameters  (7,  6,)  over  this  region.  Preliminary 
studies  (28]  of  alternative  criterion  that  consider  the  statistical  nature  of  deavage  initia¬ 
tion  and  the  volume  of  the  process  zone  support  the  use  of  a  simple,  single-point  compar¬ 
ison  for  the  present  assessment  of  size  effects. 

For  a  specified  material,  the  microscopic  failure  criterion  is  related  to  applied  7-values 
in  SENB  spedmens  through  the  procedure  illustrated  in  Fig.  2. 7-values  for  SENB  speci¬ 
mens  ( Isenb )  are  plotted  against  7-values  for  the  SSY  model  (Jssy)  th*t  impose  the  identi¬ 
cal  opening  mode  stress  at  r  -  4<3, on  the  crack  plane.  Thus,  when  Jsenb  •  Jssy  the  crack- 
tip  stresses  match  those  for  an  infinite  body,  fracture  toughness  values  obtained  in  SENB 
spedmens  under  such  conditions  are  geometry  independent.  However,  large-scale  plas- 
tidty  and  size  effects  require  significant  increases  in  Jsenb  relative  to  Jssy  to  maintain 
identical  crack-tip  stresses  as  indicated  in  Fig.  2.  When  generated  for  spcdfic  material 
properties  and  spedmen  sizes,  results  of  the  form  shown  in  Fig.  2  may  be  used  to  scale 
fracture  toughness  data  to  account  for  the  alW  ratio  by  indexing  a  given  7C  for  one  alW 
through  the  Jssy  to  obtain  Jc  for  a  different  alW.  The  applications  are  two-fold:  1 )  the  ge¬ 
ometry  dependence  of  fracture  toughness  data  for  different  a/M'ratios  may  be  eliminated 
through  correlations  to  a  single  n/M'ratio  or  to  the  SSY  model,  and  2)  commonly  available 
test  results  for  deep  notch  spedmens  ma*'  be  appropriately  scaled  for  fracture  assessments 
of  shallow  notched  structures. 


NUMERICAL  PROCEDURES 

Plane-strain  finite  element  analyses  were  conducted  to  obtain  very  detailed  resolutions  of 
the  crack-tip  stress  fields  for  the  SSY  model  and  for  SENB  spedmens.  The  numerical 
models  employ  conventional  small-strain  theory  and  address  a  wide  range  of  material 
strain  hardening.  The  remainder  of  this  section  outlines  the  details  of  the  models  and  solu¬ 
tion  procedures,  beginning  with  those  aspects  common  to  both  the  SSY  and  SENB  analy¬ 
ses. 

CONSTITUTIVE  MODEL 

The  material  model  employs  J2  deformation  plastidty  theory  (nonlinear  elastidty).  The 
uniaxial  stress-strain  curve  follows  the  Ramberg-Osgood  form  of  Eq.  (3)  with  an  addi¬ 
tional  term  linearly  dependent  on  stress,  thus 

e/e0  =  a/a0  +  a(o(o0)n  (6) 


The  total  strain  components  are  related  to  the  total  stresses  by  Eq.  (4)  with  the  addition  of 
deviatoric  and  volumetric  terms  that  vary  linearly  with  stress  according  to  Eq.  (6), 


€>j  =  ( 
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where  is  the  o'**  is  the  trace  of  the  stress  tensor  and  <5ty  is  the  Kronecker  delta.  Details  of  the 
numerical  implementation  follow  those  described  by  de  Lorenzi  and  Shih  (29]. 
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The  finite-element  analyses  for  each  configuration  are  repeated  for  three  levels  of 
strain  hardening,  high  (n  =  5),  moderate  ( n  *  10),  and  very  low  ( n  «  50).  The  parameters 
is  assigned  a  value  of  1.  Other  properties  typical  of  those  for  a  moderate  strength  structur¬ 
al  steel  are  adopted  in  the  computations,  £o  =  0.002,  o0  =  60  ksi  (414  MPa)  and  Pois¬ 
son’s  ratio,  v  =  0.3;  in  this  case  aQ  corresponds  to  the  0.2%  offset  yield  strength,  Oys  ■  Fig¬ 
ure  3  illustrates  the  relative  degrees  of  hardening  represented  by  the  selected  exponents. 

FINITE  ELEMENTS  AND  CRACK-TIP  MODELS 

Eight-noded  isoparametric  elements  with  reduced  (2x2)  Gauss  quadrature  are  employed 
in  all  of  the  models.  Reduced  integration  eliminates  locking  of  arbitrarily  shaped  elements 
once  the  incompressible  conditions  of  fully  plastic  deformation  constrain  volumetric 
changes  in  the  displacement  gradients  [30]. 

Concentric  rings  of  polygonal-shaped  elements  surround  the  crack-tips.  The  inner¬ 
most  ring  contains  elements  collapsed  into  wedge  shapes  with  side  nodes  retained  at  the 
mid-point  position.  Initially  coincident  nodes  at  the  crack  tip  are  unconstrained.  This 
model  produces  a  1/r  strain  singularity  appropriate  for  n  -*  »  and  allows  blunting  defor¬ 
mations  of  the  tip. 

Preliminary  analyses  of  all  configurations  and  hardening  exponents  were  conducted 
using  coarse  meshes  to  obtain  the  range  of  crack-dp  opening  displacements  present  at 
large-scale  yielding.  The  final  meshes  were  then  constructed  to  have  concentric  rings  of 
elements  enclosing  the  crack  tip  with  sizes  that  increase  geometrically  with  r.  The  mesh 
gradation  provides  at  least  8-10  elements  over  the  domain  26,  <r<  10d,  except  for  the 
very  earliest  stages  of  loading  when  6,  is  truly  infinitesimal. 

Stresses  are  computed  at  the  2x2  Gauss  point  locations  with  nodal  stresses  obtained  by 
a  least-squares  extrapolation  of  the  Gauss  point  values  [31].  Extrapolated  nodal  values 
are  then  arithmetically  averaged  without  regard  to  sizes  of  elements  incident  on  the 
nodes.  Comparisons  of  opening  mode  stresses,  Oyy ,  along  0  =  0  for  the  various  configura¬ 
tions  utilize  the  averaged  stresses  of  comer  nodes.  Mesh  sensitivity  studies  demonstrated 
the  adequacy  of  this  procedure. 

Crack  tip  opening  displacements  in  the  finite  element  models  are  determined  by  ap¬ 
plying  the  ±  45°  intercept  procedure  to  the  blunted  crack  tips. 

SMALL-SCALE  YIELD  (SSY)  MODELS 

Figure  4  shows  the  element  mesh  defined  for  the  circular  domain  of  radius  R  that  contains 
a  sharp  crack  with  tip  located  at  r = 0.  Symmetrical  boundary  conditions  are  applied  on  the 
crack  plane  (X  £  0,  7=0).  Loading  of  the  model  is  accomplished  by  imposing  displace¬ 
ment  increments  of  the  elastic  singular  field  for  Mode  I  on  the  outer  circular  boundary: 

A u(r,  0)  =  (1  + fcf-)i(3  -  4v- cos0 )  cos(0/ 2)  (8a) 

E  2ji 
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Av(r,  6)  - 


(8b) 
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is  increased  until  the  plastic  zone  reaches  **R/10-R/15  at  which  point  the  computed  /  is 
slightly  smaller  thanAT^l  -  v^/E. 

The  mesh  contains  2277  nodes  and  720  elements  divided  into  40  rings  of  18  equally 
sized  fans.  Element  sizes  decrease  geometrically  as  r  -*  0  to  focus  the  greatest  refinement 
near  the  crack  tip.  Analyses  with  additional  mesh  refinement  did  not  alter  the  predicted 
stresses  over  the  crack  plane. 

SENB  MODELS 

Three  a/W ratios  are  considered:  0.05, 0.15,  and  0.50.  The  finite-element  model  defined 
for  the  three-point  bend  specimen  with  a!W=  0.15  is  shown  in  Fig  5.  Conventional  plan- 
form  dimensions  are  employed,  span  S  =  4  W,  with  W=  2  in.  (5.08  cm)  used  in  the  numeri¬ 
cal  solutions.  Symmetry  conditions  permit  modeling  of  only  one-half  of  the  specimen. 
The  half-symmetric  mesh  contains  350  elements  and  1200  nodes;  very  similar  mesh  densi¬ 
ties  are  defined  for  a!W-  0.5, 0.05.  The  core  of  240  elements  for  the  crack-tip  region  is 
common  for  all  three  meshes.  The  core  elements  are  divided  into  30  rings  of  increasing 
size  with  8  equally  spaced,  A0,  elements  per  ring.  Elements  incident  on  the  crack  tip  have 
side  length  L  =  0.0002  in.  (0.005  mm).  The  much  greater  refinement  in  r  relative  to0is 
necessary  to  provide  8-10  elements  over  the  region 2<5f  <r<  10<5,  at  Iow-to-moderate 
loading  levels. 

The  mid-span  load  is  distributed  over  two  small  elements  to  minimize  effects  of  the 
local  singularity  of  an  idealized  point  load. 

SOLUTION  PROCEDURES 

Loads  are  increased  to  final  values  in  variably  sized  increments,  typically 20-40  steps,  with 
full  Newton  iterations  performed  within  each  increment  to  remove  residual  forces.  Very 
stringent  criterion  are  defined  which  ensure  convergence  of  strains  and  stresses  in  the 
third  significant  figure.  Three  to  five  iterations  are  generally  required  for  convergence 
with  the  larger  number  of  iterations  and  load  steps  needed  for  very  low  hardening  (n  =  50). 
Given  the  strain  path  independence  of  stresses  in  deformation  plasticity,  converged  values 
for  strains  and  stresses  are  invariant  of  the  load  step  sizes  used  in  the  computations.  The 
number  and  magnitude  of  the  load  steps  is  selected  to  provide  a  complete  description  of 
the  response  history  and  to  maintain  quadratic  convergence  of  the  Newton  iterations. 

/-integral  values  are  obtained  using  the  domain  integral  method  [32,33].  The  varia¬ 
tion  among /-values  computed  over  element  rings  adjacent  to  the  crack  tip  elements  and 
rings  remote  from  the  tip  is  less  than  1%  as  expected  for  the  very  refined  models  and  de¬ 
formation  plasticity. 
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Numerical  computations  were  performed  on  Apollo  workstations  using  PATRAN  [34  J 
for  mesh  generation  and  POLO- FINITE  [35]  for  analyses  and  domain  integral  computa¬ 
tions. 


COMPUTATIONAL  RESULTS 


SSY  MODELS 

The  SSY  models  provide  the  tensile  stress  distributions  on  the  crack  plane  for  comparison 
and  normalization  of  the  SENB  results.  In  nondimensional  form.  Figure  6  compares  the 
finite-element  SSY  results  for  n  =  10  with  the  HRR  values  given  by  Eqs.  (2).  SSY  results 
are  shown  for  two  levels  of  applied  loading  that  produce  plastic  zones  with  sizes 
rp/R  =  1/35  and  1/10.  For  the  larger  plastic  zone,  Kj  computed  from  the  domain  integral 
/-value  using  the  elastic  conversion  differs  less  that  2%  from  the  Kj  imposed  through  dis¬ 
placements  on  the  circular  boundaty.  The  SSY  stresses  possess  similarity  length  scales 
roo/J,  r/dt  and  reveal  a  slowly  growing  divergence  from  the  HRR  stresses.  At  r  -  4  <5,  the 
SSY  stress  for  n  =  10,  for  example,  has  decreased  to  90%  of  the  HRR  value.  Similar  plots 
constructed  for  n  —  5, 50  and  for  other  angular  orientations  display  equal  scalability  of  the 
SSY  stresses  with  ra0/J  and  r/6t . 

Normalizations  of  the  SENB  fields  require  the  SSY  stress  for  arbitrary,  nondimen¬ 
sional  positions  along  0  =  0.  For  convenience,  continuous  functions  are  constructed  to  fit 
the  finite-element  nodal  values  of  SSY  stress  which  take  the  form: 

oyy  =  o0af?ecP  (9a) 

where 

P  =  r/(J/aao€0)  (9b) 

and  a,  b,  and  c  are  curve  fitting  parameters.  This  functional  representation  is  adopted 
strictly  to  minimize  errors  in  curve  fitting;  no  attempt  is  made  to  enforce  the  HRR  field 
for  r/(Jaa0€o)  <0.0005  where  finite-element  stresses  may  be  less  reliable.  The  table  in¬ 
cluded  in  Fig.  6  provides  the  curve  fitting  parameters  for  each  strain-hardening  exponent. 

SENB  SPECIMENS 

Figures  7-11  provide  the  essential  results  of  the  finite-element  analyses  needed  to  corre¬ 
late  local  stress  fields  for  the  different^ /Bfratios.  At  increasing  levels  of  applied  load,  Fig. 
7  shows  the  variation  of  stress  component  Oyy  directly  ahead  of  crack  tip  for  a!W-  0.15  and 
n  =  10.  The  SENB  results  are  normalized  by  stresses  of  the  SSY  model  when  loaded  to  the 
same  /-value.  The  normalized  stresses  are  plotted  against  distance  ahead  of  the  crack-tip 
measured  in  multiples  of  6f .  Points  that  would  lie  within  the  blunting  zone  of  a  finite-de¬ 
formation  analysis,  r  <  2-3  6, ,  are  not  shown.  The  intensity  of  deformation  for  each  curve 
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is  indicated  through  a  reference  limit  load,  P0,  defined  for  a  perfectly  plastic  material  of 
yield  stress  ob  [936]  and  through  the  magnitude  of  CTOD  relative  to  the  crack  length. 

Figure  7  shows  that  the  SENB  stresses  steadily  decrease  below  the  SSY  values  for  in¬ 
creasing  load.  Very  similar  responses  are  obtained  for  the  other  hardening  exponents  and 
tf/IF ratios.  Agreement  between  SENB  and  SSY  stresses  increases  with  greater  degrees  of 
strain  hardening  and  larger  a/H'ratios,  i.e .,alW~  0.5  withn-  5  provides  the  closest  agree¬ 
ment  with  the  SSY  solution  over  the  complete  range  of  loading.  The  relative  "flatness"  of 
the  curves  for  P/P0  <-l  indicates  that  SENB  stresses  ahead  of  the  crack-tip  have  the 
same  spatial  variation  as  the  SSY  stresses  when  distances  are  normalized  by  r/6,  or 
r/(J/a0).  Amplitudes  of  the  SENB  stresses,  however,  are  substantially  lower  than  SSY 
values  at  the  same  applied  /. 

The  similarity  in  the  spatial  variation  of  the  SENB  and  SSY  stresses  also  applies  along 
other  rays  of  constant  6.  Figure  8  compares  nondimensional  contours  of  principal  stress 
for  the  SSY  solution  with  those  for  an  SENB  specimen  of  a  !W~  0.5  at  a  low  level  of  load¬ 
ing.  Both  the  shape  and  amplitude  of  the  contours  coincide  except  for  the  small  spike  at 
±  45°  which  is  attributed  to  the  coarser  mesh  in  the  angular  direction  for  the  SENB  mod¬ 
els  (the  element  mesh  for  the  SSY  model  has  twice  the  angular  refinement  of  the  SENB 
mesh).  The  areas  enclosed  by  the  SSY  and  SENB  contours  in  Fig.  8  agree  to  within  3%. 
Figure  9  illustrates  the  effects  of  large-scale  yielding  on  nondimensional  contours  of  prin¬ 
cipal  stress  for  n  - 10  and  a/W-  0.05, 0.5.  Although  the  contours  for  both  cases  maintain 
the  constant  shape  of  SSY,  their  size  (when  normalized  by  J)  decreases  with  increased 
loading  and  extent  of  plastic  deformation.  The  absolute  size  of  the  contour  actually  in¬ 
creases  with  J,  but  at  a  slower  rate  than  predicted  by  the  SSY  field,  even  for  the  highly 
constrained  a/W=  0.5  geometry. 

Figures  7-9  provide  two  key  results.  First,  Fig.  7  clearly  indicates  that  for  each  load 
level  a  constant  scaling  factor  is  applicable  to  correlate  the  SENB  and  SSY  stresses  at  least 
over  the  length-scale  of  4-85, .  The  scaling  factor  is  defined  by  the  ratio  of  JsenbUssy 
necessary  to  bring  the  SENB  stresses  into  agreement  with  the  SSY  stresses  over  this  re¬ 
gion.  Since  a  critical  value  ioxJssY  represents  a  size-independent  cleavage  fracture  tough¬ 
ness,  the  J senb Ussy  ratio  quantifies  the  geometry  dependence  olJc,  the  measured  fracture 
toughness  of  a  finite-size  specimen.  For  deeply  notched  specimens  with  moderate  to  high 
hardening,  JsenbUssy  approaches  unity  for  ao0/J  greater  than  »200.  For  shallower 
cracks  and  light  hardening,  this  ratio  exceeds  3-4,  i.e.,  for  a  given  Jssy,  significantly  larger 
/-values  are  needed  in  shallow  notched  SENB  specimens  to  generate  equivalent  SSY 
stresses. 

The  second  key  result  of  Figs.  7-9  is  the  insensitivity  of  the  distance  ahead  of  the  crack 
tip  and  angle  6  at  which  SENB  and  SSY  stresses  are  equated  to  define  the  JsenbUssy  ratio. 
The  simple  procedure  adopted  for  the  present  study  enforces  equality  of  stresses  at 
r  =  4<5,  on  6  =  0  to  define  the  JsenbUssy  ratio.  The  more  elaborate  scheme  which  equates 
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areas  enclosed  by  principal  stress  contours  to  define  the  Jsesb^ssy  ratio  predicts  essen¬ 
tially  identical  scaling  factors  (28). 

Figures  10-12  present  the  effects  of  a/H'and  hardening  exponent  on  the  JsenbUssy 
ratio.  Figure  2  represents  the  conceptual  form  for  these  results  and  illustrates  their  use  to 
correlate  toughness  values  between  a/W  ratios.  The  data  in  Figs.  10-12  are  presented  as 
the  effectivt  driving  force  for  cleavage,  Jssy,  plotted  against  the  apparent  (measured)  driv¬ 
ing  force,  Jsenb.  The  dashed  line  on  each  figure  represents  the  small-scale  yielding  limit, 
where  Jsenb  *=Jssy  by  definition.  The  /-values  are  normalized  with  the  remaining  liga¬ 
ment  length,  b,  rather  than  crack  length  simply  to  better  separate  the  curves.  The  flow 
stress,  rather  thano0,  is  also  used  to  normalize  /  to  make  the  results  readily  applicable  to 
other  than  power-law  hardening  materials.  The  flow  stress  for  a  Ramberg-Osgood  mate¬ 
rial  is  estimated  with  the  following  relationship: 


Ofiow  -  y  (1  + 


exp(N)] 


(10) 


where  N=*  1  In.  Equation  (10)  is  derived  by  solving  for  the  tensile  instability  point  in  Eq. 
(6),  converting  true-stress  to  engineering  stress,  and  averaging  a0  and  the  estimated  ten¬ 
sile  strength. 


Note  that  scales  on  the  figures  are  different  to  accommodate  the  full  extent  of  loading 
in  each  case.  The  curves  in  Figs.  10-12  all  agree  with  the  SSY  limit  at  low  /-values  but 
show  significant  deviation  as  /  increases.  The  divergence  from  SSY  occurs  more  rapidly 
and  at  lower  /-values  in  the  shallow  notched  specimens  and  in  low  hardening  materials. 
For  n~50  (Fig.  12),  the  effective  driving  force  saturates  at  a  constant  value,  particularly  so 
for  a/W=  0.5;  further  increases  of /do  not  affect  Jssy-  Once  a  specimen  reaches  the  satura¬ 
tion  value  of  Jssy.  stresses  ahead  of  the  crack-tip  no  longer  increase  and  the  likelihood  of 
cleavage  fracture  with  further  loading  decreases  considerably.  Cleavage  in  such  cases 
would  seem  probable  only  if  the  crack  grew  by  ductile  tearing  and  sampled  a  microstruc- 
turally  weak  feature.  The  procedure  to  correlate  toughness  values  between  a/W7 ratios  and 
to  SSY  conditions  is  graphically  illustrated  in  Fig.  2. 


APPLICATION  TO  EXPERIMENTAL  DATA 


Vr:h  the  aid  of  Figs.  10-12,  experimental  values  of  cleavage  fracture  toughness  may  be 
corrected  (or  scaled)  to  account  for  the  loss  of  constraint  due  to  large-scale  (in-plane) 
yielding.  Given  the  measured  toughness  (/c),  specimen  size,  alW  ratio,  and  material  hard¬ 
ening  characteristics,  the  equivalent  toughness  (Jssy)  for  an  infinite  size  specimen  can  be 
estimated.  Alternatively,  the  critical  /  for  one  a/W,  for  example  0.15,  can  be  estimated 
from  the  measured  Jc  for  another  a/W,  e.g.,  0.5.  Both  procedures  are  applied  here  to  the 
extensive  set  of  fracture  toughness  data  for  A36  steel  recently  reported  by  Sorem,  et  al. 
[1]- 

Sorem  tested  SENB  specimens  having  a  square  section  (B  x  B)  with  thickness  of  31.8 
mm  (1.25  in.)  for  the  A36  material.  Full  transition  curves  were  generated  for  a/W  ratios  of 
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0.15  and  0.5.  In  addition,  three-dimensional,  finite-element  analyses  were  performed  to 
support  the  conversion  of  measured  loads,  load-line  displacements  and  crack  mouth 
opening  displacements  into  applied  /-values  without  assumptions  on  the  location  of  rota¬ 
tion  points,  11-factors,  etc.  The  uniaxial  stress-strain  curve  for  the  A36  material  exhibits  a 
well-defined  yield  plateau  followed  by  relatively  high  strain  hardening.  The  power-law 
behavior  of  the  Ramberg-Osgood  model  in  Eq.  (6)  is  reasonably  followed  in  the  strain¬ 
hardening  region.  The  hardening  exponent,  n,  for  this  material  is  estimated  to  be  6. 

Figure  13  shows  Jc  test  results  for  the  A36  steel;  values  for  individual  specimens  are 
plotted  for  a/W-  0. 15  while  only  the  lower-bound  curve  is  shown  for  a!W=  0.5.  At  -43°C 
and  below,  all  but  one  of  the  a!W=  0.15  specimens  foiled  by  cleavage  without  significant 
stable  crack  extension.  A  fibrous  thumbnail  denotes  evidence  of  ductile  tearing  prior  to 
unstable  fracture.  The  shallow  notch  specimens  show  a  significant  increase  in  the  apparent 
fracture  toughness  throughout  the  transition  region.  Jc  values  for  atW-  0.5  represented  by 
the  lower-bound  curve  are  corrected  for  constraint  loss,  using  Figs.  10  and  1 1,  to  predict  a 
corresponding  lower-bound  curve  for  alW=  0.15  using  n  =  5,10  hardening  exponents. 
Figure  13  shows  that  the  lower-bound  curve  predicted  for  n  =  5  slightly  under  estimates 
the  correction  needed  for  constraint  loss  while  the  n  =  10  curve  over  corrects. 

The  typical  scatter  in  toughness  values  shown  for  the  a/H/=0.15  specimens  tested  at 
-43°C  and  -76°C  suggests  that  each  test  result  at  a  specified  temperature  should  be  cor¬ 
rected  individually  for  constraint  loss.  The  crack-tip  stress  fields  differ  significantly  at  fail¬ 
ure  due  to  the  plastic  deformation  for  the  minimum  and  maximum  toughness  values,  par¬ 
ticularly  at  -43°C.  Consequently,  the  largest  correction  for  constraint  loss  should  be  , 
applied  to  the  specimen  failing  with  the  maximum  toughness;  similarly,  the  specimen  fail¬ 
ing  with  the  lowest  toughness  requires  the  least  correction.  Figures  10  and  11  are  used 
again  to  obtain  factors  that  convert  the  measured  toughness  of  each  test  specimen  (deep 
and  shallow  notch)  to  a  constraint  corrected  SSY  value.  Linear  interpolation  between 
n-5  and  n  =  10  is  applied  to  obtain  the  factors  for  n  =  6  corresponding  to  the  A36  steel. 
Figure  14  shows  all  the  measured  toughness  values  for  these  two  temperatures  and  alW 
ratios  along  with  the  corrected  small-scale  yielding  values  (Jssy)  for  each  data  point.  The 
values  of  toughness  corrected  for  constraint  loss  at  each  temperature  agree  very  well;  dif¬ 
ferences  in  “corrected”  toughnesses  for  the  two  temperatures  reflect  the  real  differences 
in  the  toughness  property  of  the  material.  Moreover,  the  corrected  values  have  less  scatter 
than  the  uncorrected  data.  The  decreased  scatter  reduces  the  number  of  specimens  that 
must  be  tested  to  establish  specific  confidence  limits  on  lower-bound  values.  Anderson 
and  Dodds  [28]  provide  a  complete  statistical  analysis  of  the  A36  data  and  the  effects  of 
the  constraint  correction  procedure. 

CONSTRAINT  LOSS  CORRECTIONS  FOR  CTOD 

Figures  10-12  provide  the  correction  factor  for  constraint  loss  in  terms  of  the /-integral. 
Corresponding  figures  may  be  constructed  with  CTOD  ( 6t )  rather  than  /  as  the  fracture 
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characterizing  parameter.  The  HRR  relationship  between  J  and  6,  defined  in  Eq.  (5)  is 
generalized  using  the  material  flow  stress  into  the  form 


J  =  moflowd,  (11) 

where  the  coefficient  m  is  often  termed  a  constraint  factor  and  is  dependent  upon  the 
specimen  geometry  and  material  flow  properties.  Table  1  provides  m-values  for  the  three 
a/FKratios  and  hardening  exponents  considered  in  the  present  work.  Flow  stress  values  for 
the  power-law  hardening  model  are  again  defined  by  Eq.  (10).  Linear  regression  is 
applied  to  plots  of  Jvs.  Oftowdi  to  obtain  m  for  each  case.  The  geometry  independent  m- 
values  for  the  SSY  model  also  shown  in  Thble  1  are  obtained  by  simple  conversion  of  the 
HRR  dn  values,  which  are  found  here  to  hold  for  the  SSY  models.  The  factors  in  Table  1 
may  be  applied  to  express  the  axes  of  Figs.  10-12  in  terms  of  dt  rather  than  J. 


a/W 

Strain-Hardening  Exponent 

n  =  5 

n  =  10 

n  =  50 

0.05 

1.78 

1.50 

1.30 

0.15 

1.85 

1.57 

1.45 

0.50 

2.04 

1.79 

1.69 

SSY 

2.10 

1.71 

1.37 

Table  1.  Constants 'm  ’  to  Convert  J  to  CTOD,  <5,  =  //(monow) ,  for  Ramberg-Osgood 
Material  Models 


DISCUSSION 

Finite-element  analyses  for  the  SENB  specimens  convincingly  demonstrate  the  effects  of 
in-plane  plastic  deformation  on  the  amplitude  of  crack-tip  stresses.  The  loss  of  constraint 
(stress  triaxiality)  caused  by  impinging  of  global  plastic  deformation  on  the  crack  tip  low¬ 
ers  the  local  stresses  in  the  SENB  specimens  relative  to  small-scale  yield  values,  when  the 
SENB  and  SSY  fields  are  compared  at  the  same  J.  The  divergence  of  SENB  and  SSY  stress 
fields  is  accelerated  by  decreasing  the  a/JFratio  and  lowering  the  degree  of  strain-harden¬ 
ing.  For  a/W=  0.05  and  0. 15,  agreement  with  the  SSY  stress  field  is  achieved  only  at  very 
low  levels  of  loading  and  vanishingly  small  crack-tip  plastic  zones.  Consequently,  cleav¬ 
age  fracture  toughnesses  obtained  from  shallow  notch  specimens  in  the  transition  region 
will  exhibit  a  strong  geometry  dependence.  Limits  of  the  type  in  Eq.  (1)  to  assure  size  inde¬ 
pendent  Jc  values  for  shallow  notches  would  be  so  severe  as  to  be  impractical.  Further- 
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more,  the  analyses  show  that  SENB  specimens  of  very  low  hardening  materials  containing 
deep  notches  \a!W=  0.5)  also  exhibit  some  loss  of  constraint  that  influences  the  local 
stress  field  once  large-scale  plasticity  is  achieved,  i.e.,  deformations  imposed  following 
development  of  a  plastic  hinge. 

Principal  stress  contours  near  the  crack  tip  for  the  SENB  specimens  are  self-similar 
and  closely  match  the  shape  of  contours  for  the  SSY  solution,  when  distances  from  the 
crack  tip  are  normalized  using  ra0/J  or  r/6, ,  i.e.,  the  similarity  length  scales  for  SSY.  Am¬ 
plitudes  of  the  SENB  stresses,  however,  scale  with  /in  a  geometry  dependent  way  which  is 
substantially  different  from  that  predicted  by  the  geometry  independent  SSY  model. 
These  observations  suggest  a  crack-tip  field  for  the  SENB  specimens  which  separates  the 
principal  stresses  into  multiplicative  functions  in  the  same  spirit  as  HRR,  such  that 

o?ENB{r,  6)  =  hi(a/W;  aa0/J ;  n)  •  a,SSY  i  -  1, 2. 3  (12) 

where  the  full  SSY  fields  have  the  form 


a,jSY(r,  6)  =  a0  *  a,j(ro0/J-,  n\&)  (13) 

Both  hi  and  dy  are  dimensionless  functions  of  their  arguments.  In  the  crack-tip  region 
where  the  plastic  component  of  strain  greatly  exceeds  the  elastic  component,  the  /i2,j  func¬ 
tions  are  related  through  the  von  Mises  equivalent  stress  to  hj,  which  can  be  determined 
from  a  finite-element  analysis.  The  most  relevant  aspect  of  Eqs.  (12,13)  is  the  separation 
of  SENB  stresses  into  a  geometry  and  material  dependent  amplitude  (or  constraint)  func¬ 
tion,  hi,  which  does  not  vary  with  r  and  0 .  When  the  loading  is  sufficiently  small  that  ampli¬ 
tudes  of  crack-tip  fields  for  the  SENB  specimens  approach  those  of  SSY  conditions. 
Moreover,  off* -*  off** as  r-*0.  Although  closed-form  expressions  fordo  in  Eq. 
(13)  are  not  yet  available,  resolution  appears  feasible  through  additional  analyses  of  the 
SSY  model  for  a  broader  range  of  material  parameters.  The  framework  of  Eqs.  (12)  and 
(13)  represents  a  logical  extension  of  Parks  [20]  recently  proposed  concept  of  organizing 
elastic-plastic  crack-tip  fields  into  “similar”  classes  of  increasingly  complex  form  which 
degenerate  to  HRR  under  appropriate  restrictions. 

The  existence  of  crack-tip  fields  that  follow  the  form  of  Eqs.  (12-13)  permits  a  rela¬ 
tively  straightforward  construction  of  the  relationship  between  SENB  and  SSY  stresses 
for  a  given  material  and  geometry  (a,  W)  as 

aSENB^aSSY.  HiMm)  (14) 

JSSY 

Finite-element  analyses  of  the  type  presented  in  Section  4  for  the  SENB  specimens  pro¬ 
vide  estimates  for  the  H,  functions  (the  Hj  functions  defined  using  stresses  on  the  symme¬ 
try  plane,  for  example,  are  represented  by  the  curves  plotted  in  Figs.  10-12). 

The  SENB  stress  fields  presented  in  Section  4  as  determined  by  the  finite-element 
analyses  are  for  idealized  plane-strain  conditions.  When  the  specimen  thickness  is  finite. 


14 


however,  through-thickness  constraint  can  be  considerably  less  than  plane-strain  condi¬ 
tions.  Three-dimensional  analyses  of  SENB  specimens  with  in-plane  mesh  refinement 
comparable  to  that  used  for  the  2-D  models  in  the  present  study  are  very  rare  due  to  the 
enormous  computational  expense.  One  such  analysis,  recently  performed  by  Narishimhan 
and  Rosakis  [37],  provides  some  insight  regarding  the  thickness  required  to  maintain 
nearly  plane-strain  conditions  through  the  thickness.  They  analyzed  a  SENB  specimen 
having alW  -  0.4  where  the  crack  length  and  remaining  ligament  length,  b,  were  three  and 
six  times  the  thickness,  B,  respectively;  thus  thickness  was  the  governing  dimension.  Very 
low  hardening,  n  =  22  was  specified  for  their  analysis.  Figure  15,  constructed  from  the  Na¬ 
rishimhan  and  Rosakis  results,  shows  the  stress  normal  to  the  crack  plane,  relative  to  the 
mid-thickness  value.  Data  for  three  load  steps  are  plotted  which  correspond  to  Boa/J  ra¬ 
tios  of  235, 103,  and  26.3.  The  relative  distance  ahead  of  the  crack  tip  for  each  case  lies  in 
the  range  of  2-4  <5,.  At  the  lowest  7-value,  the  opening  mode  stress  is  nearly  constant 
through  the  thickness  except  very  close  to  the  free  surface.  At  the  intermediate  load  step, 
the  stress  remains  relatively  constant  over  the  middle  40%  of  the  thickness;  while  at  the 
largest  7-value,  the  stress  decreases  continuously  from  the  mid-thickness  value  toward  the 
free  surface.  Narishimhan  and  Rosakis  also  performed  a  2-D  plane-strain  analysis  of  the 
SENB.  The  3-D,  crack-tip  stresses  at  mid-thickness  displayed  very  good  agreement  with 
the  2-D,  plane-strain  values  for  the  two  lower  7-values  shown  in  Fig.  15.  Thus,  the  deeply 
notched  SENB  specimen  seems  to  maintain  nearly  plane-strain  constraint  through  a  sig¬ 
nificant  portion  of  the  thickness  for  Baa/J  ratios  at  least  into  the  150-200  range. 


By  combining  the  plane-strain  analyses  of  the  present  study  with  the  3-D  analyses  of 
Narishimhan  and  Rosakis,  specimen  size  limits  to  assure  geometry  independent  cleavage 
fracture  can  be  established,  at  least  for  deeply  notched  specimens.  Figure  16  shows 
Jsenb/Jssy  ratios  for  the  a/W=0.5  configuration  as  determined  by  the  finite-element 
analyses  of  Section  4.  The  curves  become  relatively  flat  and  approach  unity  when 
the  aoaow/J  ratio  exceeds  approximately  200,  although  the  rate  at  which  each  curve  ap¬ 
proaches  the  SSY  limit  depends  on  the  hardening  exponent.  Thus,  the  recommended  size 
requirements  for  deeply  notched  specimens  are: 


B,  b,  a 


2007c 

Oflow 


(15a) 


or 


B,b,a  >  300<5C 


(15b) 


These  requirements,  which  would  also  apply  to  deeply  notched  compact  specimens,  guar¬ 
antee  fracture  toughness  results  that  are  nearly  size  independent,  but  only  when  cleavage 
occurs  without  significant  stable  crack  growth.  Eq.  (15a)  is  eight  times  more  severe  than 
the  size  requirements  in  E8 13-87  but  not  nearly  as  severe  as  E399-83.  Consider,  for  ex¬ 
ample,  a  material  with  7C  =  200  kPa  m,  oys  =  450  MPa,  and  =  500  MPa.  The  minimum 
thickness  required  for  a  valid  A/c  test  is  570  mm  (22.4  in),  while  a  10  mm  (0.39  in)  thick  ' 
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specimen  satisfies  the  E8 13-87  requirement.  An  80  mm  (3.15  in)  thick  specimen  is  re¬ 
quired  to  satisfy  Eq.  (15a). 

While  the  present  work  has  focused  on  effects  of  the  in-plane  geometry  on  fracture 
toughness;  the  effect  of  thickness  requires  considerably  more  study,  especially  the  interac¬ 
tion  of  small  tf/W ratios  and  thickness.  Figure  15  shows  that  deeply  notched  test  specimens 
can  maintain  nearly  plane-strain  conditions  at  mid-thickness  to  relatively  high  7-values, 
but  the  size  of  the  plane-strain  region  decreases  with  plasticity.  It  should  be  possible  to 
define  an  effective  thickness,  which  equals  the  actual  thickness  for  small-scale  yielding, 
but  which  decreases  with  J  for  large-scale  yielding. 

The  effect  of  prior  ductile  crack  growth  also  requires  further  study.  The  results  pres¬ 
ented  here  apply  strictly  to  stationary  cracks.  Ductile  crack  growth  affects  cleavage  tough¬ 
ness  in  at  least  two  ways:  (1)  the  stress  fields  ahead  of  a  growing  crack  differ  from  those  of  a 
stationary  crack;  and  (2)  the  growing  crack  samples  more  material  which  increases  the 
likelihood  of  encountering  a  microstructural  feature  that  triggers  cleavage.  Nevertheless, 
it  seems  reasonable,  at  least  for  a  first-approximation,  to  correct  cleavage  fracture  values 
for  constraint  loss  using  the  results  of  Figs.  10- 12  provided  the  amount  of  ductile  growth  is 
less  than  the  size  of  the  finitely  deformed  blunting  zone,  i.e.,  A  a  <  *2  <5, . 

SUMMARY  AND  CONCLUSIONS 

Very  detailed  finite-element  analyses  for  shallow  and  deeply  notched  SENB  specimens 
have  supported  the  development  of  a  framework  to  quantify  the  effects  of  constraint  loss 
due  to  in-plane  plastic  deformation  on  cleavage  fracture  toughness.  The  numerical  results 
provide  the  relationship  between /and  crack-tip  stress  fields  when  the  specimen  geome¬ 
try  and  extent  of  plasticity  invalidate  the  asymptotic  HRR  and  small-scale  yielding  (SSY) 
solutions.  Spatial  variations  of  crack-tip  stresses  for  the  SENB  specimens  closely  match 
those  for  the  small-scale  yielding  model  of  crack-tip  behavior,  when  distances  from  the 
crack  tip  are  normalized  as  ra0/J.  However,  amplitudes  of  the  SENB  stress  fields  vary 
with  J  in  a  substantially  different  manner  than  for  small-scale  yielding.  Similar  to  the  form 
of  the  HRR  solutions,  these  results  suggest  a  separation  of  the  crack-tip  stresses  into  an 
amplitude  (or  constraint)  function  which  multiplies  the  SSY  fields,  where  the  material  and 
geometry  dependent  amplitude  term  is  determined  through  finite-element  analyses  of  the 
type  performed  in  this  study.  Consequently,  the  relationship  between  SENB  and  SSY 
crack-tip  stresses  can  be  expressed  simply  through  functions  of  the  Jsenb/Jssy  ratios  for  a 
given  material.  Graphical  forms  of  these  functions  are  derived  from  the  finite-element 
solutions  and  when  applied  to  an  A36  steel  are  shown  to  eliminate  size  effects  on  cleavage 
toughness  data  in  the  transition  region.  Moreover,  the  large  scatter  in  toughness  values 
inherent  in  experimental  testing  is  significantly  diminished  through  the  constraint  correc¬ 
tion  process. 

The  numerical  results  also  demonstrate  that  crack-tip  stresses  for  deeply  notched 
SENB  specimens  with  some  strain-hardening  closely  approach  small-scale  yielding  con- 
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ditions.  Fracture  toughness  data  obtained  under  such  conditions  should  be  independent 
of  specimen  geometry.  To  assure  J  controlled  deavage  fracture  m  deeply  notched  bend 
specimens,  the  results  of  this  study  indicate  that  all  relevant  dimensions  of  a  specimen 
should  exceed  200  JcIOfom-  This  requirement  is  eight  times  more  strict  that  current  stan¬ 
dards  for  Ju  testing  but  substantially  less  stria  than  requirements  for  Aj,  testing 
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Fig  1.  Schematic  /-temperature  transition  curve  for  shallow  and  deeply  notched  SENB  specimens.  Typical  extents  of  plastic  zones  are 
shown  for  shallow  and  deeply  notched  specimens. 


SSY 
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Fig  2.  Schematic  of  relationship  between  the  /-integral  values  in  SENB  specimens  for  deep  and  shallow  cracks  with  those  of  SSY 
conditions  which  generate  equivalent  opening  mode  stresses  ahead  of  the  crack  tip. 


K|  Displacement  Field 
Imposed  on  Boundary 


21 


Finite-element  mesh  for  the  small-scale  yielding  computations. 


Typical  finite-element  mesh  for  plane  strain  analyses  of  the  SENB  specimens. 


Comparison  of  normalized  crack  opening  stress  in  plane  strain  small-scale  yielding  vs.  HRR  solution  for  hardening 
exponent  n  =  10  at  two  levels  of  loading.  Thble  of  curve  fitting  parameters  for  small-scale  yielding  response  for  hard¬ 
ening  exponents  n  -  5,10,  and  50. 


10,  alW =0.5 
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Fig  8.  Comparison  of  principal  stress  contours  in  small-scale  yielding  and  in  a  SENB  specimen  for  a!W- 0.5,  harden 

ing  exponent  n  - 10. 


0.020 


Fig.  10.  Relationship  between  the  normalized  7-integral  values  in  SENB  specimens  with  those  of  SSY  conditions  which  generate  equiva¬ 

lent  opening  mode  stresses  ahead  of  the  crack  tip  for  hardening  exponent  n  -  5. 
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Fig  11. 
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Temperature  (°C) 

Fig  13.  Experimentally  determined  critical  /-values  for  an  A36  steel  showing  the  estimated  lower-bound  curves  for  a!W=  0.15  gener 

ated  from  the  a/W=  0.5  curve  using  the  constraint  loss  correction  procedure. 


A36  STEEL 

=  31.8  mm  {1.25  in. 
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Fig  14.  Comparison  of  experimentally  determined  critical  7-values  for  an  A36  steel  with  data  values  corrected  for  constraint  loss. 


Narishimhan  &  Rosakis) 


Comparison  of  crack  opening  stress  through  the  thickness  of  a  deeply  notched  SENB  specimen  relative  to 
centerplane  values. 
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Fig  16.  Ratio  of  /-values  in  deeply  notched,  SENB  specimens  to  /-values  in  small-scale  yielding  that  generate  the  same  crack  opening 
stresses  for  increased  levels  of  applied  loading  and  hardening  exponents  n  -  5,  10,  and  50.  Proposed  size  requirement  for  Jc 
testing  to  assure  geometry  independent  toughness  measures  and  current  size  requirement  for  //c  are  indicated. 
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